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Abstract. Inspired by the multiverse scenario, we study a heterotic dark energy model in 
which there are two parts, the first being the cosmological constant and the second being the 
holographic dark energy, thus this model is named the AHDE model. By studying the AHDE 
model theoretically, we find that the parameters d and £lhde are divided into a few domains 
in which the fate of the universe is quite different. We investigate dynamical behaviors of 
this model, and especially the future evolution of the universe. We perform fitting analysis 
on the cosmological parameters in the AHDE model by using the recent observational data. 
We find the model yields Xmin = 426.27 when constrained by Planck + SNLS3 + BAO + HST, 
comparable to the results of the HDE model (428.20) and the concordant ACDM model 
(431.35). At 68.3% CL, we obtain —0.07 < Dao < 0.68 and correspondingly 0.04 < QhdeO < 
0.79, implying at present there is considerable degeneracy between the holographic dark 
energy and cosmological constant components in the AHDE model. 
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1 Introduction 

The holographic dark energy model (HDE) [1, 2] was motivated by the holographic principle 
[3], as one of promising models to solve the nature of dark energy [4, 5]. The basic idea behind 
the HDE is that our universe is a sense finite and can be described by a two dimensional 
spherical holographic screen, thus there must be finite size effects, and one of these effects is 
the contribution to the zero point energy, depending on the size of the screen. Parametrically, 
this contribution assumes the form 


Phde = 3d 2 M pl L 2 , (1.1) 

where d is a dimensionless parameter to be determined by experiments, M p i is the reduced 
Planck mass, and L denotes the size of holographic screen. For convenience, we work in the 
natural units, where H = c = 1. In [2], one of the present authors suggested to choose the 
future event horizon of the universe as the size of the holographic screen, given by 

r+°° dt 

Rh = a -. (1.2) 

Jt a 

This choice not only gives a reasonable value for the dark energy density, but also leads to 
an accelerated expansion. 

The holographic dark energy (HDE) model based on Eq. (1.1) and Eq. (1.2) has proven 
to be a promising dark energy candidate. In the original paper [2], Li showed that the HDE 
can explain the coincidence problem. In [6], it is proven that the model is perturbatively 
stable. Other studies show that the model is in good agreement with the current cosmological 
observations [7]. Thus, the HDE model becomes one of the most competitive and popular 
dark energy candidates, and attracts a lot of interests [8]. 
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It remains quite a mystery that to date all the papers on the HDE assume that dark 
energy is dominated by the HDE given by Eq.(l.l). In retrospect, this can be explained 
only by the reason that all the authors believed that the universe dominated by the HDE is 
the unique universe thus the form of dark energy is also unique. In the multiverse scenario, 
however, our observable universe is only one out of numerous universes, and the cosmological 
constant is one of the physics constants varying from one universe to another. Thus, it is not 
reasonable to simply assume it be vanishing. 

Of cause, the modern multivese scenario was motivated by the problem of dark energy. 
One of the present authors (ML) has recently been converted into a believer of the multiverse 
scenario by a quite different problem, namely the Fermi paradox. He now believes that the 
correct answer to the Fermi paradox is that our human being is the only intelligent being in 
our galaxy, possibly the only intelligent being in our universe, since it is really very difficult 
for an intelligent being to appear, the only reason for us to appear is that our universe is one 
out of numerous universe and it just happens that our universe is lucky enough, this is an 
anthropic answer to the Fermi paradox. 

If for whatever reason that the multiverse scenario is true, then it is natural that the 
cosmological constant indeed is a constant to be determined by observations. On the other 
hand, the HDE on a general ground must be present too, according to the holographic prin¬ 
ciple. Thus, dark energy must consist of two parts, the first is a constant, the second is of the 
finite size effect. Indeed, in a calculation of the photon contribution to the zero point energy 
[9], it is found that in addition to the usual UV divergent part, there is a second divergent 
part proportional to L~ 2 where L is the radius of the de Sitter space. The usual quartically 
divergent part can be absorbed into the cosmological constant, the second part is the same 
form of the HDE. Thus, in general 


Pir — A + bMp t L 2 + b\L ^ + ... . (1-3) 

It just happens that in our universe, the first term and the second term are comparable, and 
the third term is much smaller and can be neglected completely. 

In this paper, we shall study this heterotic model of dark energy, in which there are the 
cosmological constant and the HDE. Thus we have 

Pde PA T Phde■ (1-4) 

We shall name this model AHDE model. This model raises interesting theoretically questions: 
Is the CC positive or negative? If it is positive, whether it is greater or smaller than the HDE? 
Is the future of our universe dominated by the CC or by HDE? Will the big rip happen or 
not? In the rest of this paper, we shall try to answer these questions. 

This paper is organized as follows. We study the AHDE model theoretically in Sect.II, 
and find that the parameters d and Qhde are divided into a few domains in which the fate 
of the universe is quite different. In Sect.Ill, we fit the model to the combined Planck + 
SNLS3 + BAO + HST and Planck + SNLS3 + BAO + HST + SDSS - Lya datasets, and 
present the fittings results in Sec. IV. Many interesting issues, including the ratio of HDE 
and the equation of state (EoS) of the AHDE are discussed. Some concluding remarks are 
given in Sec. V. In this work, we assume today’s scale factor ao = 1, so the redshift z satisfies 
z = 1/a — 1. The subscript “0” indicates the present value of the corresponding quantity 
unless otherwise specified. 
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2 AHDE Model: Theoretical Analysis 


In this section, we will write down the basic equations for the AHDE model in a non-flat 
universe and study the fate of the universe in this model. 

2.1 HDE with Cosmological Constant 

In a spatially non-flat Friedmann-Robertson-Walker universe, the Friedmann equation can be 
written as 

( 2 . 1 ) 


^MpiH 2 — Pdm + Pb + Pr + Pk + Pdei 


where p k = —3is the effective energy density of the curvature component. In the AHDE 
model, the dark energy density is 

Pde = PA + Phde = M^A + 3 M^id 2 R h 2 . (2-2) 

For convenience, we define the fractional energy densities of the various components, i.e., 


o _ ^ r> _ Pde O _ Phde 

— yj2 2 ’ ^ — 5 ^ ^hde — 5 

i± CL p c p c 


Q _ PA Q _ Pdm /-v _ Pb Q _ Pr 

— 5 dm — ? ^6 — 5 — > 


Pc Pc Pc pc 

(2-3) 

where p c = 3M^H 2 is the critical density of the universe. The subscripts, “/c”, “he”, “hde”, “A”, 
“dm”, “6” and “r” represent curvature, total dark energy, holographic dark energy, cosmological 
constant, dark matter, baryon and radiation, respectively. By definition, we have 

hde + ^A + Hrfm + + Q k = 1. (2-4) 

The energy conservation equations for the components in the universe take the forms 

Phde T 3H(phde T Phde ) — 0, (2-5) 

PA = o, Pdm + 3Hp dm = 0, p b + 3Hp b = 0, p r + 4Hp r = 0, p k + 2Hp k = 0. (2.6) 

Combining Eq. (2.5) and Eq. (2.6) together, we can obtain the form of phde, 


2 H 


1 1 


Phde — g 2 P c P c 3 A + 3 dfc + PA ■ 


(2.7) 


Substituting phde into Eq. (2.5), follow the similar procedure in Ref. [10], we get a differential 
equation of H and Qhde'- 


2(14/ide — 1)^ A &hde + H(3£lhde ~' 3 + 3J4a + ^k — ^ r ) — 0. 

From the energy density of the HDE in Eq. (1.1), we have 

L — d 


Hy/^lhde 

Following Ref. [11], in a non-flat universe, the IR cut-off length scale L takes the form 

L = ar(t), 


( 2 . 8 ) 


(2.9) 
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( 2 . 10 ) 







and r(t) satisfies 


r r{t) 


dr 


r+oo 

Jt 


dt 


Jo \/l — kr 2 Jt a {t) 

By carrying out the integration, we find 


r{t) = 


1 


VW\ 


sinn 


1 


+°° dd\ 

t a) ^/\k\ 


sinn 


where the function sinn(x) is dehned as 

{ sin(x) k > 0; 
x k = 0; 

sinh(x) k < 0. 

Equation (2.10) leads to another equation about r(t), namely, 

, , L d 

r(t) = = 


a V^hdeHa 


Combining Eqs. (2.12) and (2.13) yields 


. y/W\d 

\k / — = arcsmn- 

lt a 



"s/^hde&H 

Taking derivative of Eq. (2.14) with respect to t, one gets 

^ 4 hde 


Z^hde + + H 


VhdeH 2 k 


d 2 


2.2 Evolution Equations of E{z) and Qhdei z ) 


+°° da \ 

a(t) Htf)' 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


Combining Eq. (2.8) with Eq. (2.15), we eventually obtain the following two equations gov¬ 
erning the dynamical evolution of the AHDE model in a non-flat universe, 


1 dE(z) _ f ihde ( 3Qa + flfc - Q r — 3 1 

E{z ) dz 1 + 2 l 2 Q hde 2 


^ 4 hde 

~dT 


+ ^A 


(2.16) 


d^hde “^Qhdei 1 ^ hde ) 


^4 hde q 1 314 a T ^4 r 

~cp~ k 2 “ 2(1 - n hde ) 


(2.17) 


dz 1 + 2 

where E(z) = H(z)/Ho is the dimensionless Hubble expansion rate. Notice that we have 

o — ^fco(l + z ) 2 n _ S4r°(l + z) 4 n _ S4ao n _ S4&o(l + z) 3 fr > 

^k\Z) — t~i/ \o j il r yZ) — . 2 ? ^ 6 (‘ 2 ') — t-i/ \o 5 ( 2 . 18 ) 


E(z) 


E{z) 


E{zf 


E(z) 


and the fractional density of dark matter is given by 14 dm(z) = 1 — £ 4 a :(~) — ^hde{z) — £ 4 a ( z ) ~ 
£l r (z) — kib(z). Equations (2.16) and (2.17) can be solved numerically and will be used in the 
data analysis procedure. 
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2.3 Dark Energy Equation of State 

The EoS of the HDE takes the form [12] 


^hde 


Phde 


1 2 
3 ~~ 3 


^hde 


+ 


Phde 3 3 V cP 

So according to the partial pressure law, the EoS of the total dark energy is 

^ _ Pde _ Phde T PA _ ^hde^hde ~ ^A 

Pde P\ide T PA ^hde T 

Obviously the property of w^ e is closely related to values of d and 


(2.19) 


( 2 . 20 ) 


2.4 The Fate of the Universe in AHDE Model 

For convenience, we transform Eq. (2.1) into the following form 

(1 ~ $lhde)H~ = tlmoHya 3 + Q, r oHQa 4 + a ^^AoHq, (2.21) 


and define 

/(o) = fl m o-£fo« 1 + QroHyCL 2 + + QaoHq a~. 

Let x = log a, then we obtain 


H = 


/(«) 


a 2 (l - sihde) 


( 2 . 22 ) 


(2.23) 


d , Qfide 

— in -—- 

dx 1 \ihde 


+ s m|/ ( a)| 


2 

d 


hde (1 Qhde ) 


kd 2 

/(a)' 


(2.24) 


Eq. (2.24) tells us how the HDE evolves with a. For simplicity, in this section we only 
study the k = 0 case. Then Eq. (2.24) becomes 


d , ^hde 

— m -—— 

dx 1 ^Ifide 


+ s M/ ( a)| 



(2.25) 


This equation can not be solved exactly. But for the purpose of studying the fate of the 
universe, we can introduce a good approximation. 

During a period in which /(a) is dominated by a single term on the right hand side of 
Eq. (2.22), we can use a constant k r to approximate hr |/(o)|. Then the Eq. (2.25) becomes 


d , hde 

— in -—— 

dx 1 D/ide 


2 

d 




hde 


(2.26) 


It is easy to see that the cosmological constant term on the right hand side of Eq. (2.22) 
will dominate when scale factor a evolves to infinity. Thus we get k r = 2 in this limit. 
Fortunately, it can be proved that a will always evolve to infinity. The proof is given in the 
appendix. 
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y/Clh 




(a) 0 < y/0.hde 0 < 1 (b) 0 < y/0.hdeO < d (c) d < y/Qhde 0 < 1 



V Qhde 


V Qhde 


y/^hde 

- • - < 


- • - < 


- • - 

- • - 


d 1 Id Old 


(d) 1 < %/fi/ideO 


(s) d <1 \/it h.rU.O 


(f) 1 < V^hdeO < d 


Figure 1. The evolution of y/Cl^de with different initial conditions. 


With k r = 2, we can rewrite the Eq. (2.26) and get its solution as follows: 


dClhde 

dx 

X + X\ 


cjiV^hde di)(l Clhde'jClhdei 


d 


2d -2 
d 


In 


1-y/n 


hde 


i 


0 d 2 - 1 ) 


In 


\/Clhde d 


In Clhde T 


+ 


2d+ 2 


In (1 + 

hde ) i 


(2.27) 


(2.28) 


where x\ is the constant of integration. Then let us exhibit how the HDE evolves: 

The first three pictures in Fig.l show how y/Clhde evolves with the initial condition 
y/ Clhde o < 1- In both picture(la) and picture(lb), y/ Clhde o < d, the right hand of Eq.(2.27) 
is negative, so ^Me 0 will decrease with x. Thus the right hand of Eq.(2.27) will always 
be negative. So y/ClhdeO will keep decreasing to 0. According to Eq.(2.28) when y/ClhdeO —> 
0 + ,x —>• +oo. Thus once y/Clhde 0 < 1, either d < 1 or d > 1, according to Eq.(2.23), H will 
always tend to a constant, which means space-time will be de Sitter in the future. In picture 
(lc), d < y/Tlhde o < 1, the right hand of Eq.(2.27) is positive, so y/ Clhde o will increase with 
x. Thus the right hand of Eq.(2.27) will always be positive, y/Cl hde 0 will keep increasing to 
1, according to Eq.(2.28) when^/Qft.cfeo —t C~,x —> +oo. So this solution describes a HDE 
dominated universe with big rip. 

The last three pictures in Fig. 1 show how y/Clhde evolves with the initial condition 
V Clh.deO > 1- In picture (Id), y/ Clhde o > d, the right hand of Eq.(2.27) is negative, so when 
d < 1, the y/Cl hde o will keep decreasing with x all the way up to 1. According to Eq.(2.28) 
wheni/H/j^eo —> l + ,x —> +oo. So picture(ld) describes a HDE dominated universe with 
big rip. In picture(le) (or in pic.ture(lf)), d < yj Clhde o (or d > y/Clhde o) the right hand of 
Eq.(2.27) will be negative (positive), so the y/ClhdeO will keep decreasing (or increasing) to 
d. According to Eq.(2.28), when^/fl/i^eo — > d, x —> +oo. So they both describe an universe 
whose space-time will be de Sitter in the future. 
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3 The Observational Data and Methodology 

We explore cosmological constraints on the AHDE model with the most recent observational 
data. For comparison, we will also present the fitting results of the original HDE and ACDM 
models. Data used in our analysis include: 

• The SNLS3 combined sample [13, 14], consisting of 472 SNIa, combining the results of 
two light-curve fitting codes SiFTO [15] and SALT2 [16]. 1 We follow the procedure of 
[25] and perform a complete analysis of the systematic errors. The y 2 function takes 
the form 

XSNLS 3 = A ^ T ■ C_1 ' ( 3 - 1 ) 

where C is a 472 X 472 covariance matrix capturing the statistic and systematic un¬ 
certainties, and Arri = rr^ — n^ mo d is a vector of model residuals of the SNIa sample, 
with tub the rest-frame peak B band magnitude of the SNIa and m mot i the predicted 
magnitude of the SNIa, given by 


mmod = 51og 10 T> L - a(s - 1) + /3C + M, (3.2) 

where T>l is the Hubble-constant free luminosity distance, the stretch s is a measure of 
the shape of SN light-curve, C is the color measure for the SN, and a, /3 are two nuisance 
parameters characterizing the stretch-luminosity and color-luminosity relationships, re¬ 
spectively. Following [13], we treat a and /3 as free parameters of y 2 function. 

• The Planck “distance priors”provided in [26], which are extracted from Planck first year 

[27, 28] observations. The data include the baryon component ujb = flbh 2 , the “acoustic 
scale” l a = 7rr(z*)/r s (z*), and the “shift parameter” R = r(z Jf ), where z* is 

the redshift to the photon-decoupling surface [29], r(£*) is our comoving distance to £*, 
and r s (z*) is the comoving sound horizon at z*. The distance priors provide an efficient 
summary of the CMB data in regards to dark energy constraints [30]. 

• The BAO data including the measurement of r s /Dy at z = 0.106 from 6dFGS [31], 
the isotropic measurement of Dy/rd at z = 0.32 from the BOSS DR11 LOWZ sample 
[32], the anisotropic measurement of D\/r L \ and Hrd at z = 0.57 from the BOSS DR11 
CMASS sample [32], and the improved measurements of Dy/r s at z = 0.44, 0.60, 0.73 
from the WiggleZ Dark Energy Survey [33]. Here rj is the comoving sound horizon at 
the “drag” epoch when the baryons are “released” from the drag of the photons [34], and 
Dy is a volume averaged distance indicator similar to the angular diameter distance 
Da [35], 

• The Hubble constant measurement Hq = 73.8 ± 2.4km/s/Mpc from the WFC3 on the 
HST (Hubble Space Telescope) [36]. 

• The high-redshift BAO measurement from the Quasar-Lya-forest cross-correlation of 
the BOSS DR11 of SDSS-1II [37], namely ajj ,7 a^ 3 = 1.025 ± 0.021; ay and aj_ are 

1 It should be mentioned that, previous studies on the SNLS3 data sets [17] found strong evidence for the 

redshift-dependence of color-luminosity parameter /3, and this conclusion has significant effects on parameter 

estimation of various cosmological models [18-21]. In addition, different light-curve fitters of SNIa can also 
affect the results of cosmology-fits [22-24]. But in this work, for simplicity, we just adopt the most mainstream 
recipe of processing SNLS3 data and do not consider the factors of time-varying (3 and different light-curve 
fitters. 
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Table 1 . Fitting results for the AHDE model 



Planck + SNLS3 + BAO + HST 

Planck + SNLS3 + BAO + HST + SDSS - Lya 

Parameter 

Best fit 

68.3% limits 

Best fit 

68.3% limits 

H 0 

d 

^A0 

0.1405 

73.3 

0.003 

0.61 

-0.0012 

0 141 q+ uuul::a 

U.JU1D_ 0 0025 

70.5l{j 

0 f;7n+ 0 - 32 0 

U.O/U_ 0 .180 

0 07+ 0 - 61 
u ’ u ‘ -0.14 

0 nnq7+ 0 0037 
u - uuo ' -0.0059 

0.1375 

71.6 

0.001 

0.63 

-0.0063 

0 1400 +U ' UU2,S 

u.muu_ 0 0025 

70.5t{i 

0-308 

0.42lgf 0 
—0.0013l°;°°§« 

Y 2 • 

Amm 

426.27 


431.79 



defined as a\\ = > a -L = [D^ly/rlfd ’ where 2 : = 2.34, the fiducial values 

[-D//(z)/rd] fid and [Da(z)/ rd] &d are 8.708 and 11.59 respectively. 

In the following context, we will use “SNLS3”, “Planck”, “BAO”, “HST” and “SDSS — Lya” 
to represent these five datasets. 

We combine the above data sets to perform y 2 analyses. Since SNLS3, Planck, BAO, 
HST and SDSS-Lya are effectively independent measurements, the total y 2 function is just 
the sum of all individual y 2 functions: 


2 _ 2 ,2 ,2 ,2 , 2 

Xtotal — XSNLS3 + Xpianck T XBAO + XHST + XSDSS-lycr 


(3.3) 


In our work, for a detailed investigation, we do fittings with two datasets: Planck + SNLS3 + 
BAO + HST and Planck + SNLS3 + BAO + HST + SDSS — Lyct, respectively. 

The AHDE model has two dark energy parameters d, and Including four other 
cosmological parameters Q m h 2 , ojb , and h. and two nuisance parameters a, (3 characterizing 
the systematic errors of the SNLS3 dataset [13], the full set of free parameters in our analysis 
is 

P — {H m /i , iObi h : c , D/j, oi , /?}. (3.4) 

In this work, we numerically solve Eq.(2.16) and Eq.(2.17) to obtain background evolutions 
of the AHDE model. The values of D r o, for simplicity, are determined from the 7-yr WMAP 
observations [38], 

Q r0 = D 7 q( 1 + 0.2271W//), ^ 7 o = 2.469 x HT 5 /?" 2 , N eff = 3.046, (3.5) 


where 7 represents photons, and N e ff is the effective number of neutrino species. 

We modify the public available CosmoMC package [39] to explore the parameter space 
using the Markov Chain Monte Carlo (MCMC) algorithm. All the parameters listed in Eq. 
(3.4) are fitted simultaneously. 


4 Dynamical Behaviors and the Cosmic Expansion History 
4.1 Fitting Results 

In Table 1, we give best-fit parameters as well as 68.3% confidence limits for constrained 
parameters. The results show that the spatial curvatures are close to zero in both cases (the 
68.3% confidence limits are \Qko\ < 0.007 and |I < 0.004, respectively). Thus the results 
are impressively consistent with a spatially flat universe. Table 1 also gives the constraint on 
the parameter 
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Figure 2. Marginalized likelihood distribution of fl m h 2 , fl A o, d constrained by Planck + SNLS3 + 
BAO + HST (solid lines) and Planck + SNLS3 + BAO + HST + SDSS — Lya (dashed lines) datasets. 





d 


Figure 3. Marginalized 68.3% and 95.4% CL contours of VL\Q-Q, m h 2 , d-fl m h 2 and flAo-d planes 
constrained by Planck + SNLS3 + BAO + HST (solid lines) and Planck + SNLS3 + BAO + HST + 
SDSS — Lya (dotted lines) datasets. 


Planck + SNLS3 + BAO + HST: -0.07 < O A0 < 0.68 (68.3% CL); 

Planck + SNLS3 + BAO + HST + SDSS - Lya: 0.32 < O A0 < 0.67 (68.3% CL). 

Combine with the constraint results of Q m h 2 and Hq listed in Table 1 together, we get 
the corresponding constraint on Qhde o : 

Planck + SNLS3 + BAO + HST: 0.04 < O Me0 < 0.79 (68.3% CL); 

Planck + SNLS3 + BAO + HST + SDSS - Lya: 0.06 < il hdeQ < 0.41 (68.3% CL). 

Fig. 2 shows the marginalized likelihood distributions of parameter d, £l m h 2 and fi A o 
constrained by Planck+SNLS3+BAO+HST and Planck+SNLS3+BAO+HST+SDSS—Lya 
datasets, respectively. The 68.3% and 95.4% contours of H A o-H m h 2 , d-Q m h 2 and H A o -d planes 
are plotted in Fig. 3. 

The results show that, compared with the Planck + SNLS3 + BAO + HST dataset, the 
Planck + SNLS3 + BAO + HST + SDSS — Lya dataset makes a more tighter constraints on 
d and H A o parameters. We also find that a smaller value of D m h 2 and a bigger value of H A o 
is favored by the Planck + SNLS3 + BAO + HST + SDSS — Lya dataset. 
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Table 2. The Xmin s i AAICs and ABICs of the ACDM, HDE and AHDE models, obtained by using 
the Planck + SNLS3 + BAO + HST and Planck + SNLS3 + BAO + HST + SDSS — Lya datasets, 
respectively. 


Model 

Planck + SNLS3 + BAO + HST 

Planck + SNLS3 + BAO + HST + SDSS - Lya 

Amin 

AAIC 

ABIC 

a. mm 

AAIC 

ABIC 

ACDM 

431.35 

0 

0 

438.22 

0 

0 

HDE 

428.20 

-1.15 

3.03 

438.19 

1.97 

6.15 

AHDE 

426.27 

-1.06 

7.27 

431.79 

-2.43 

5.92 


There is a degeneracy between d and Oaoj which can be seen in Fig. 3. The reason is 
that both the cosmological constant and HDE are good candidate for explaining the feature of 
cosmic acceleration revealed by current observational data. Therefore, when we combine the 
HDE and cosmological constant components together (thus AHDE model), we may probably 
get the degeneracy. 

In addition to the cosmological consequence of the AHDE model, we are also interested in 
its comparison with the ACDM and the original HDE models. Therefore we also perform the 
X 2 analysis of the ACDM and the HDE models by using the same datasets. To assess different 
models, here we adopt the Akaike information criteria (AIC) [40] and Bayesian information 
criteria (BIC) [41], defined as 

AIC = Xmin + 2fc, BIC = Xmira + ^ In-W) (4.1) 


where k is the number of free parameters, and N is the number of data points used in the 
fits. A model with smaller AIC (BIC) is more favored. 

Table 2 shows the Xmin s > AICs and BICs of the ACDM, HDE and AHDE models. 
Notice that the values of the AIC and BIC themselves are not interesting, thus we only list 
the difference between the AHDE (HDE) and ACDM models, i.e., 


AAICmodei — AIC mo dei AICacdm> ^BIC moc iei — BIC mo( iei BICacdm (4.2) 


Compared with the ACDM and HDE models, the AHDE model provides a better fit to the 
data. For Planck + SNLS3 + BAO + HST dataset, the AHDE model reduces the Xmin s by 
amount of 5.08 (1.93) compared with the ACDM (HDE) model. While for Planck + SNLS3 + 
BAO + HST + SDSS — Lya dataset, the AHDE model reduces the Xmin s by amount of about 
6.4 compared with the other two models. By adopting the AIC and BIC, Table 2 also shows 
that, compared with ACDM model, the AHDE model is slightly favored by AIC. But both the 
AHDE and HDE models are not favored by BIC, though these models have slightly smaller 
Xmin s than the ACDM model. 


4.2 The Expansion History 

It is worth investigating the cosmic expansion history of the AHDE model by the fitting 
results. 
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Figure 4. Reconstructed evolution history of H(z) (95.4% CL) in AHDE model, constrained by the 
Planck + SNLS3 + BAO + HST (red solid lines) and Planck + SNLS3 + BAO + HST + SDSS — Lya 
(black dashed lines) datasets respectively. 

In Fig. 4, we plot the reconstructed evolution history of H(z) (95.4% CL) in AHDE 
model, constrained by the Planck+SNLS3+BAO+HST and Planck+SNLS3+BAO+HST + 
SDSS — Lya datasets, respectively. We find that, in low redshift region, the reconstructed 
evolution history H(z) of the Planck + SNLS3 + BAO + HST and Planck + SNLS3 + BAO + 
HST + SDSS — Lya datasets are almost the same. However, in the high redshift region, the 
constraint of the Planck + SNLS3 + BAO + HST + SDSS — Lyct dataset is much more tighter 
than that of the Planck + SNLS3 + BAO + HST dataset. It is clear that this feature is due 
mainly to the SDSS-Lya data at redshift z = 2.34. As revealed by [42, 43], the high redshift 
datasets play a big part in constraining the cosmic expansion history. 



Figure 5. Reconstructed evolution history of H(z) (95.4% CL) constrained by Planck + SNLS3 + 
BAO + HST + SDSS — Lya dataset in AHDE model(black solid lines) . The corresponding results for 
HDE model is also plotted (red dashed lines) for comparison. 


- 11 - 




For a comparison, it is also of value to compare the AHDE model with the original HDE 
model. Fig. 5 shows the reconstructed evolution history of H(z) (95.4% CL) in the AHDE 
and the original HDE model constrained by Planck + SNLS3 + BAO + HST + SDSS — Lya 
dataset. We find that the reconstructed H(z) of AHDE and HDE models have negligible 
difference in the low redshift region. However, in the high redshift region, the H(z) in AHDE 
model has slightly lower value, which should be due mainly to the existence of a cosmological 
constant component in the model. 

4.3 Equation of State 

In this subsection we discuss the EoS w, which is believed to be the most important marker 
of the properties of dark energy. 

Fig. 6 shows the reconstructed evolution history of w(z ) at 0 < z < 2.5(68.3% and 
95.4% CL) constrained by Planck + SNLS3 + BAO + HST + SDSS — Lya dataset. It shows 
that w slightly cross -1 from above roughly at the current epoch. However, in the past we 
have w slightly bigger than -1, which can be viewed as a feature of diluted holographic dark 
energy. This behavior is consistent with the results shown by Fig. 5. 

As mentioned above, the dynamical evolution of dark energy have not be confirmed by 
the current observational data. Our results is consistent with this statement. 

From Fig. 4, we can conclude that, if we want to break the degeneracy between the 
HDE and cosmological constant components, one way is to get more observational data at 
high redshifts. 



Figure 6. Reconstructed evolution history of w(z) at 68.3% (black solid lines) and 95.4% (red dashed 
lines) CL in AHDE model, constrained by the Planck + SNLS3 + BAO + HST + SDSS — Lya dataset. 

5 Conclusion 

In this work, we study the AHDE model in which there are two parts, the first being the 
cosmological constant and the second being the holographic dark energy. This model has 
similar dynamical equations with the original HDE model, except a cosmological term. By 
studying the AHDE model theoretically, we find that the parameters d and Q,hde are divided 
into a few domains in which the fate of the universe is quite different. 
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Using the Planck + BAO + SNLS3 + HST and Planck + BAO + SNLS3 + SDSS — Lya 
datasets, we investigate the dynamical properties and cosmic expansion history of the AHDE 
model. The results shows that the goodness-of-fit of the AHDE model are Xmin = 426.27 
(Planck+SNLS3+BAO+HST) and Xmin= 431 - 79 (Planck+SNLS3+BAO+HST+SDSS-Lya) 
which is smaller then the results of the original HDE model (Planck + SNLS3 + BAO + 
HST:428.20; Planck + SNLS3 + BAO + HST + SDSS — Lya:438.19) and the concordant 
ACDM model (Planck + BAO + SNLS3 + HST:431.35; Planck + SNLS3 + BAO + HST + 
SDSS — Lya:438.22) obtained using the same datasets. Especially when constrained by the 
Planck+SNLS3 + BAO + HST + SDSS — Lya dataset, The Xmm °f AHDE model shrinks more 
than 6, compared with both the HDE and AHDE model. Thus, the AHDE model provides a 
nice fit to the cosmological data. 

For parameter fl A o, the 68.3% confidence level constrained by the Planck + SNLS3 + 
BAO + HST and the Planck + SNLS3 + BAO + HST + SDSS - Lya dataset is -0.07 < P A0 < 
0.68 and 0.32 < U A o < 0.67, respectively. This gives the corresponding components of the 
holographical dark energy, namely 

Planck + BAO + SNLS3 + HST: 0.04 < n hde0 < 0.79; 

Planck + SNLS3 + BAO + HST + SDSS - Lya: 0.06 < tt hdeQ < 0.41. 

We also find that there is degeneracy between the cosmological constant and the holo¬ 
graphic dark energy component when constrained by current cosmological observations. By 
reconstructing the evolution of the EoS of dark energy, we find the AHDE mainly differs from 
the original HDE model at high redshift (as shown in Fig. 5). 

From the constraint results by Planck + SNLS3 + BAO + HST + SDSS — Lya dataset, it 
shows that if we want to break the degeneracy between the HDE and cosmological constant 
components, one way is to get more observational data at high redshifts. 
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6 Appendix: Proof for log a not having a maximum in AHDE model 

In this appendix, we prove that x does not have a maximum as a function of time. Since 
x does not have a maximum in either ACDM or HDE model, here we only consider AHDE 
model, which includes the coexistence case of HDE and cosmological constant. 

In a flat universe, the Friedmann equation is 

3 MpiH 2 = p dm + p d + p r + p A + Phdet (6-1) 
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Let us note that 


9 (®) = ^; ln !/(«)!, 


then g(x) = 


2Q A0 e 2x - n m0 e~ x - 2Q r0 e 

^Aoe 2x + fl m oe~ x + H r0 e 


-lx 


-lx 


According to Eq.(2.25), we get 

d^lhde 


^ [ j \Z^hde 5(^)]^/ide(l ^ hde )• 


( 6 . 2 ) 

(6.3) 


(6.4) 


We first consider the cases that Hao > 0 or Hao + H m o + H r o <0. In these cases, g{x) 
is a bounded function in the region (xo,+oo). If x has a maximum x m . H would approach 
zero when x goes to x m . According to Eq.(6.1), £lhde would approach infinity. However, 
Eq.(6.4) shows that, once y/Vl hde > %g(x) and Q hde > 1, ^V^hde < 0. Thus Q hde would 
never approach infinity. So x has no maximum in these cases. 

Another case is Hao < 0 and Hao + Sl m o + H r o > 0. In this case 


3x c > 0 : O, A0 e 2Xc + Q m0 e Xc + fi r0 e 2Xc = 0. 


(6.5) 


So Eq.(6.4) has a singularity at x = x c . Similar to the previous analysis, x does not have 
a maximum in the region (x'o, x c ). Now we prove that x c is not a maximum of x. If it is, 
g(x) would approach negative infinity and &hde would approach infinity when x goes to x c . 
However, once £lhde > 1, d L's/^hde < 0. Thus Qhde would never approach infinity. So x c is 
not a maximum of x. When x > x c . we can see that Ha + H m + H r < 0. This is just the case 
we have discussed. So x also has not a maximum when x > x c - 

In summary, we have proved that x does not have a maximum. So we can use a constant 
k r to approximate ^ In \f(a)\ when we study the fate of the universe. 
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